^ 



> 

en 

00 






IPMU-1 3-0078 

UT-13-14 

April 5, 2013 



Recent Result of the AMS-02 Experiment and 

Decaying Gravitino Dark Matter in Gauge Mediation 

en 

o 

(N 

O^ Masahiro Ibe*^"'^\ Sho Iwamoto*^^\ Shigeki Matsumoto*^^\ 

■^ Takeo Moroi*^'^'^) and Norimi Yokozaki*^'') 

Oh 

Q-i ( \ 

dj ^"^'ICRR, University of Tokyo, Kashiwa, Chiba 277-8568, Japan 



^^^Kavli IPMU, University of Tokyo, Kashiwa, Chiba 277-8568, Japan 
^'^'^ Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 



(^ Abstract 

cn 



The AMS-02 collaboration has recently reported an excess of cosmic-ray 
positron fractions, which is consistent with previous results at PAMELA and 
Fermi-LAT experiments. The result indicates the existence of new physics 
phenomena to provide the origin of the energetic cosmic-ray positron. We 
pursue the possibility that the enhancement of the positron fraction is due to 
the decay of gravitino dark matter. We discuss that such a scenario viably fits 
into the models in which the soft SUSY breaking parameters are dominantly 
from gauge-mediation mechanism with superparticle masses of around 10 TeV. 
Our scenario is compatible with 126 GeV Higgs boson, negative searches for 
SUSY particles, and non-observation of anomalous FCNC processes. We also 
point out that the scenario will be tested in near future by measuring the 
electric dipole moment of the electron and the lepton flavor violating decay of 
the muon. 



1 Introduction 

The AMS-02 collaboration has recently released their first result of the cosmic-ray 
positron fraction [1]. The anomalous excess over the expectation based on the sim- 
ple comic-ray propagation models has been seen again, and its energy spectrum is 
consistent with the previous results of the PAMELA [2] and Fermi-LAT [3] experi- 
ments. The AMS-02 collaboration has also reported a fact that the positron flux at 
the energy region of the excess shows no appreciable anisotropy to date. This fact 
may indicate that the positrons are not due to some astrophysical activities on the 
galactic plane but are from exciting dark matter phenomena in the halo, though it 
is clearly premature to make a definite statement because of limited statistics. 

When the excess of the positron fraction is interpreted as a dark matter signal, 
a decaying gravitino with its mass being 0{1) TeV could be a promising candidate 
for dark matter [H \5\ (for earlier discussions on decaying gravitino dark matter see 
references j6]). The lifetime of the gravitino is then required to be (9(10^^) sec, 
which is realized by introducing a tiny i?-parity violation. The introduction of this 
violation is in fact favored from the viewpoint of cosmology; we can evade the serious 
gravitino problem [7j because the required violation is large enough to allow other 
sparticles to decay before the era of Big-Bang Nucleosynthesis (BBN). On the other 
hand, the violation is small enough not to wash out the baryon asymmetry of the 
universe created in the early universe, so that models with the decaying gravitino 
are consistent with the successful leptogenesis scenario |8j. 

In this letter, we revisit the model of decaying gravitino dark matter in light of 
the recent results of the AMS-02 experiment and the discovery of the Higgs particle 
at the LHC experiments [U [10]. As we show, the desirable gravitino mass and the 
observed Higgs boson mass at around 126 GeV can be explained simultaneously in 
the models of direct gauge mediation with the messenger masses just below the scale 
of the Grand Unified Theory (GUT). We also show that the gravitino mass at around 
0{1) TeV allows us to construct a simple model of the i?-parity breaking which is 
appropriate for the favored gravitino lifetime, 0(10^^) sec. 

This letter is organized as follows. In section[2} we reassure that the decaying 
gravitino with a mass at around 1 TeV well fits the anomalous excess of the positron 
fraction observed at the AMS-02. In section|3j we discuss the models with gauge me- 
diation at the GUT scale which explain both the favored gravitino mass of (9(1) TeV 
and the observed Higgs boson mass, rrih — 126 GeV. We also discuss the flavor chang- 
ing neutral current and CP-violation in the lepton sector. In section]!], we give a 
simple model of i?-parity breaking which leads to the appropriate gravitino lifetime. 
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The final section is devoted to our conclusions. 

2 Positron Fraction 

We first discuss the excess of positron fractions reported by the AMS-02 collaboration 
and its implication to the decaying gravitino dark matter. There are several ways 
to introduce R-parity violating interactions which make the gravitino meta-stable. 
In this letter, we consider the simplest possibility, namely the violation through the 
LHu operator as an example, where L and H^ are the lepton doublet and the up- 
type higgs superfields, respectively. The relevant part of the superpotential and soft 
SUSY breaking terms are therefore given by as follows; 

w = ^lHuHd + ^i\HuU, (1) 

Aoft = {B^i HuHd + B[^\HuLi + h.c.) - L\m\ijLj - mjjjHd]'^, 

where H^ is the down-type higgs superfield. Then, the gravitino dark matter decays 
into a Z boson plus a neutrino, a Higgs boson plus a neutrino, and a W boson plus 
a charged lepton with the relative ratio of about 1:1:2, as was explicitly shown in 
reference [5]. 

We next summarize our procedure to calculate the positron fraction R = $e+/($e+-|- 
$e-); where <l>e+ and $e- are positron and electron fluxes, respectively. The fluxes 
consist of the contribution from the decaying gravitino dark matter and the back- 
ground contribution: $e+ = ["^e+JoM + ["^e+Jbkg- For the contribution from the dark 
matter, we have solved the diffusion equation in order to take account of the effect 
of electron/positron propagations inside our galaxy. The energy spectrum of the 
electron/positron from the dark matter /e± evolves as [TT] 

^-^'^^=K{E) [VVe±(i?,r)] + ^[6(E)/,±(E,f)]+Q(E,f), (2) 

where the function K is expressed as i^ = i^o -^Gev with -EceV being the energy in 
units of GeV and 6 = 1.0 x 1{)~^^ Eq^^ [GeV/sec]. In our numerical calculation, 
we have fixed the parameters 5 and Kq to be 0.55 and 5.95 x 10^'^ [kpc^/Myr], 
respectively. The diffusion zone is assumed to be a cylinder with a radius of 20 [kpc] 
and a half- height of 1 [kpc], and we set /e± = at the boundary. The fluxes from 
the dark matter are then obtained as [$e±]DM(-£') = {c/ An) fe±{E ^Tq) . 

The source term Q expressing the production of primary electrons/positrons from 
decaying gravitino dark matters is given by the following formula. 



QiE^r)- P^^^' 



dN, 



"^3/2 7-3/2 [ dE 

2 



(3) 

decay 



Positron Fraction 



B 

O 
Q. 

e 



m3^2=l'0TeV 
T3^2=10^^sec 



^elec;bkg=0-035E 
^posi;bkg=0-0030E 



■3.24 




E [GeV] 



Figure 1: The positron fraction expected in the decaying gravitino dark matter scenario 
with the biUnear term of the R-parity violation. The soUd hne corresponds to the L2HU 
interaction, while the dashed line to the LiH^- As a reference, the fraction without dark 
matter contribution is drawn with a dotted line. The recent results from the AMS-02 
experiment is shown with the red data points. 



where T3/2 and 7713/2 are the hfetime and the mass of the gravitino, respectively. In 
addition, pdm is the dark matter mass density of our galaxy for which we adopt the 
NFW profile pDuir) = p^ir-Q / r){rc + TqY / {tc + r)"^ [12], where p© ^ 0.4 [GeV/cm^] 
is the local halo density, Tc — 20 [kpc] is the core radius, and r© ~ 8.5 [kpc] is the 
distance between the galactic center and our solar system. In the above expressions, 
[dNf,±/dE]^cca.y is the energy distributions of the electron/positron from single decay 
process, and are calculated by using PYTHIA6 package p!3]. 

The background contributions are approximated to the form [$e±]bkg = C± ■ 
{Eq^yY^, where [$e±]bkg is in the unit of [GeV cm^ sec sr]"-*^. Parameters for elec- 
trons are determined with fitting against the result of total electron flux reported 
by PAMELA collaboration [H] to be (C_,(5__) = (0.035,-3.24). For positrons, 
C+ = 0.088C_ and 5+ = 5^ - 0.26 are used. 

The result is shown in Fig. [1} As a benchmark point, the mass and lifetime of 
the dark matter are fixed as (wi3/2,T3/2) = (1.0 TeV, 10^^ sec). Two patterns of the 
gravitino decay are shown: the solid line (the dashed line) shows the case where the 
gravitino exclusively decays into the second (first) generation leptons, i.e., R-parity 
is broken only in the second (first) generation leptons. 



First it should be noted that the data points below 10 GeV are not taken into 
consideration because the fraction in the region is subjected to the solar modula- 
tion [I5]. Also, fractions around 10 GeV could be modified by the change of the 
background contribution. Considering these two uncertainties, the AMS-02 results 
are explained well by the solid line, i.e., when the gravitino decays into the second 
generation. On the contrary, when the decay is into the first generation leptons, the 
slope becomes steeper due to the electrons produced primarily, and it is difficult to 
explain the AMS-02 result. 

Decays of the gravitino through the LH^ operator, on the other hand, also pro- 
duce other particles such as anti-protons and gamma-rays. Furthermore, high-energy 
electrons from the decay of the gravitino produce gamma-rays through the inverse 
Compton scattering with cosmic microwave backgrounds and infra-photons from star 
lights [T6] . Since no excesses are reported so far at observations of these particles, 
models involving the decaying gravitino are severely constrained. First, let us con- 
sider the constraint from the observation of cosmic-ray anti-protons. According to 
reference [TTj, with use of the most conservative limit (dark matter decay only), it 
turns out that the lifetime of the gravitino should be longer than 10^^ sec. when the 
gravitino mass is around 1 TeV. If we used shallower profiles than the NFW one, the 
limit is expected to be milder [18]. Next, we consider the constraint from gamma-ray 
observations. According to the most conservative limit (dark matter decay only) in 
reference [19], the lifetime has to be longer than about 10^^ sec. In any case, the 
lower limit on the lifetime is almost the same as that required to explain the excess 
of positron fractions. An attractive solution to this problem is the use of the LLE 
operator instead of LHu-, where we do not have to worry about the constraint from 
the observation of cosmic-ray anti-protons, while that from gamma-ray observations 
becomes much milder than the case of LH^ |20j . 

Finally, let us comment on the mass of the gravitino. The decaying gravitino 
should be heavier than 520 GeV in order to explain the rightmost data point of the 
AMS-02 result, which is for the energy of 260-350 GeV. Meanwhile, if the gravitino 
mass is heavier than 1 TeV, the lifetime which can explain the result tends to be 
shorter than 10^^ sec. since the number density of the gravitino is smaller. This is 
somewhat disfavored for the constraints from the gamma-ray and anti-proton fiuxes. 
The AMS-02 result indicates that the mass of the decaying gravitino dark matter 
should be ~ 1 TeV, which provides precious information to model building behind 
the physics of the dark matter. 



3 Gravitino Dark Matter in Gauge Mediation 

In the previous section, we have seen that the AMS-02 result is well explained in 
the model with unstable gravitino dark matter with the mass of ~ 1 TeV. However, 
such a model is usually severely constrained by the flavor violating processes, like the 
K-K mixing, lepton flavor violating (LFV) processes, and so on. In particular, if the 
gravity mediation contribution dominates the soft SUSY breaking sfermion masses, 
the off-diagonal elements of the sfermion mass-squared matrices are expected to be 
unsuppressed compared to the diagonal ones. If so, our scenario is not compatible 
with the constraints from the flavor violating processes. In addition, in the light of 
the recently observed Higgs mass, the scenario gives too small Higgs mass (unless 
the tri-linear scalar coupling is relatively large); the stop masses should be as large 
as ~ 10 TeV to realize the lightest Higgs mass as large as 126 GeV [21]. 

If the gauge-mediation contribution gives an extra contribution to the soft SUSY 
breaking parameters, the above problems may be solved. In particular, if the gauge 
mediation contributions to the sfermion masses are significantly larger than the grav- 
itino mass and are as large as ~ 10 TeV, we obtain viable scenario as we see below. 

We parameterize the messenger sector as 

"mess I, -'"mess ~r -'^mess" J^i^i; \^) 

where \l/j('l'j) are messenger fields in fundamental (anti- fundamental) representation 
of SU{b) GUT gauge group. The indices i runs to 1 — A^mess, with A'mess being the 
number of the messengers. The gravitino mass is bounded from below: 

p 

^ ^ mess /r\ 

The equality is satisfied in the direct gauge mediation models; in such a model, the 
messenger scale is around the GUT scale taking msoft ~ 10 TeV. 
The bino, wino and gluino masses are given by 

(6) 
(7) 



Mb 


r-^ 


N "^ 

-^ ' mess A 

An 


( F 

^ mess 


\-'"mess 


Mw 


r-^ 


N "' 

-^ ' mess . 

An 


/ F 

-* mess 


\ -'^mess 


M, 


r-^ 


N "' 

-' ' mess . 

An 


/ F 

^ mess 


\ -^►^mess 



where ai, a2 and aj, are gauge coupling constant squareds of 5'f/(3)c, SU{2)i and 
f/(l)y divided by An, respectively. (Here, we take a normalization of U{1)y as the 



SU{5) GUT normalization.) In addition, the messenger-scale values of the sfermion 
masses are 
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where Q,lf,D are squarks and L,E are sleptons. The soft masses of the up- and 
down- type Higgses {Hu and Ha) are same as m|. 

The mass spectrum we anticipated can be easily obtained. For instance, taking 
^mess = 10^^ GeV, -Fmess/^mess = 3000 TcV, and A^eff = 1, the low-energy values 
of gluino and squark masses are abut 20 TeV while the slepton masses are ~ 10 
TeV. Furthermore, assuming the direct gauge mediation, the gravitino mass becomes 
about 1 TeV. (Notice that the gravitino mass as large as 1 TeV can be realized with 
smaller value of Mmess if we do not consider the direct gauge mediation.) 

Let us estimate the rates of flavor and CP violations which are important probes 
of the SUSY particles even if the masses of superparticles are quite large [221 123112^. 
First, we consider the leptonic ones on which the experimental bounds are expected 
to be drastically improved. Even though all the superparticles are relatively heavy, 
the rates of LFV processes may become sizable. This is because the off-diagonal el- 
ements of the sfermion mass-squared matrices are expected to be of (^(mg/g), which 
are ~ 1 % of the diagonal elements in the present setup. If the masses of slep- 
tons, Higgsinos, and the electroweak gauginos are 10 TeV while the 1-2 elements of 
the slepton mass-squared matrices are 10^^ of the diagonal elements, for example, 
Br{fi -)■ 67) is estimated to be 1 x 10"^^, 6 x 10~^^, and 2 x 10~^^, for tan/3 = 5, 10, 
and 20, respectively. Thus the rate of the /i — )■ 67 process can be below the current 
bound. In addition, using the fact that Br{^ — )■ 67) is approximately proportional 
to tan^ /3, Br{^ — )■ 67) may be within the reach of MEG-upgrade experiment (which 
is expected to reach Br{fi — )■ 67) ~ 6 x 10~^"^ [25]) in particular when tan/3 is large. 

Another important check point is the so-called ex parameter, which often gives 
the most stringent constraint on supersymmetric flavor and CP violations. Assuming 
that off-diagonal elements of the squark mass-squared matrices are given by (1 TeV)^, 



for example, the SUSY contribution to the ex parameter becomes smaller than the 
present bound (1 x 10^'^ [221 EZj) if the mass scale of the colored superparticles is 
larger than ~ 20 TeV. (Here, the phases in the MSSM parameters are tuned to 
maximize the SUSY contribution to ex-) Such a value of the colored superparticle 
masses can be easily realized as we have mentioned. 

We also comment on the electric dipole moment (EDM) of the electron. Assuming 
that the supergravity contribution provides the B parameter of 0(7713/2), we expect 
that the CP violating phase in the /x parameter (in the bases where the VEVs of the 
Higgs bosons and the gaugino masses are real) is sizable. Taking the masses of non- 
colored superparticles to be 10 TeV, for example, the electron EDM is estimated to 
be 0.3 X 10~^^, 0.6 x 10"^^, and 1 x 10"^'' ecm for tan/3 = 5, 10, and 20, respectively, 
which are around the current bound [27jj. (Here, we have tuned the phase of the 
H parameter to maximize the electron EDM.) The future experiments, which may 
improve the bound by ~ 3 orders of magnitude |2Sj , have a good chance to observe 
the electron EDM. 

Now, we estimate the size of /i' and soft masses required for the gravitino decay 
with ~ 10^® sec. In the case that R-parity is broken by bi-linear operator, we need a 
small VEV of a left-handed sneutrino, inducing chargino, neutralino/lepton mixing. 
Then, the gravitino decays into the SM gauge bosons and leptons through the mixing. 
The sneutrino get a VEV trough the potential 

£soft 9 BfL,Hu-m%Hl + h.c.-ml^^^\L\^, (14) 

where Bf = —tiB^ + tkB'j^jjidik and mf = {5m\)iktl — eiSmfj^. The rotation param- 
eters ej are determined by the ratio of /i and /i'j as e^ = yUj/yU, and rn^^g^ = m\ are 
the soft masses induced by gauge mediation. Other mass parameters, S, B\ ^J\5mW 
and \ \5m\j | of ~ 1 TeV are induced by the supergravity effects. The sneutrino 
VEV is given by 

Ki = —^ ^{Bf sin /3 - m • cos /3) /mcMSB • (15) 

For the gravitino of ~ 1 TeV, the life-time of 10^^ sec. is explained by Kj ~ 10~^° [5]. 
This requires e^ to be ~ 10^^, i.e., the coefficient of the R-parity violating bi-linear 
term /i^ should be ~ 10^^ GeV {jj, ~ 10 TeV in the parameter space of our interest). 
This /x'j ~ 10^^ GeV can be explained consistently with the seesaw mechanism; /x' 
is induced by the small vacuum expectation value of the right-handed sneutrino 
through the operator LHu{Nr)- An explicit model realizing the small {Nr^ is 
discussed below. Note that the NLSP also decays into SM particles through the 



mixing with the hfe-time much less than 1 sec, and hence, the BBN constraint can 
be avoided [7]. 

4 A Model of R-parity Violation 

In the previous section, we have shown that the small value of the R-parity violating 
parameter of /x' = O(10~^) GeV is needed to explain the AMS-02 result. Here, we 
show that such a value of the /i' parameter may arise in the model with a discrete 
R-symmetry, as we show below. 

Let us consider a model with a discrete i?-symmetry, %^b. (see Table. [I]). With 
the given charge assignments, the relevant superpotential terms are given by, 

W = VuH^ 10 10 + yd,eHdlO 5* + y^H^ 5* Nr + fiH^H^ + i/A/0iV^/2 

''-<P'NR + ^<P' + yxx(<p' + CH^H^H,-cJ^] , (16) 



Mk Kl V "MpL " "^ ^^MpL 



where and X are gauge singlets and y's and c's are dimensionless constants,*^ We 



assume here that there is a spontaneous discrete i?-symmetry breaking sector which 



generates 1713/2 = (Wo) /Mp^ 7^ 0,*^ It should be noted that the above superpotential 
breaks U{1)b~^l and i?-parity explicitly, although the usual i?-parity violating terms 
LLE, LQD, JJUD and LH^ are suppressed by 'L^r. 

From the above superpotential, we find a stable vacuum at around, 

(0) = {^w^^\ = {cwm,i,M^^fl\ (17) 
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For cw = 0{1), dlvM = C(l) m3/2 = C(l) TeV, we find 

(0) = 0(10^^) GeV , 
(Nr) = C(10-^)GeV . (19) 

'^^ Here, we absorbed the allowed mass term of 0^ by the shift of X. We also assume that the 

size of /i is controlled by another symmetry than Z57J symmetry. 

'^^ For example, the discrete i?-symmetry is spontaneously broken in the supersymmetric Yang- 
Mills theories by the gaugino condensations 29 . If X appears in the gauge kinetic functions of the 
Yang-Mills theories, we need several gaugino condensations to stabilize X. 



Table 1: The charge assignments under the discrete Z5/J symmetry are presented here. We 
have used SU{5) GUT representations for the MSSM matter fields, i.e. 10 = {Ql, Ur, Eh), 
5* = {Dji, Li) and Nr the right-handed neutrino. 
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At around this vacuum, we obtain the masses of the right-handed neutrino, 

Mn oc (0) = C(10^^) GeV , (20) 

which are appropriate for the see-saw mechanism [30]. Furthermore, we also obtain 
the effective bi-linear R-parity violating term. 



/^ 
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0(10"'^) GeV 



(21) 



which is appropriate for the decaying gravitino dark matter discussed in the previous 
section. 



In addition to the superpotential terms in Eq. (16), there are i?-parity breaking 
terms such as, 

1 



W 



Mk 



<j){Wo)H^5* 



(22) 



which gives a similar contribution to the LHu term from the (iV/j) LH^ term. The 
trilinear i?-parity breaking terms have, on the other hand, the charge —1 under Z5R, 
and hence, they are proportional to {'m3/2/Mpi^)^^'^ = (9(10~^^). Thus, the trilinear 



i?-parity breaking terms generated by the vacuum expectation values in Eq. (17) are 
negligibly small. 



5 Summary 

In this letter, we have shown that the anomalous increase of the positron fraction 
recently observed by the AMS-02 experiment can be well explained in the unstable 
gravitino scenario if 77^3/2 ^ 520 GeV and T3/2 — 10^^ sec. Even with such a value of 
the gravitino mass, we argued that the dangerous SUSY contribution to the flavor 
violating processes may be suppressed enough if the gauge mediation mechanism 
provides the dominant contribution to the sfermion masses of ~ 10 TeV. We have 
also proposed a model to give the R-parity violating interaction required to explain 
the AMS-02 signal. 
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Here, we concentrated on the case where the gravitino decays into lepton and 



electroweak gauge boson (or Higgs) with the interaction given in Eq. (14). With 
other types of R-parity violating interaction, however, gravitino may decay into 
different final states. In particular, if the LLE-tjpe R-parity violating operator 
exists in the superpotential, the gravitino decays into purely leptonic final state, 
with which the constraints from the anti-proton flux can be easily avoided. More 
detailed analysis on the related issues will be given elsewhere. 
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